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' Abstract 



We present a fully backreacted D3-D7 supergravity solution dual to the Klebanov-Strassler cascad- 
5^ \ ing gauge theory coupled to a large number of massive dynamical flavors in the Veneziano limit. 

The mass of the flavors can be larger or smaller than the dynamically generated scale. The solution 
is always regular at the origin of the radial coordinate and as such it can be suitably employed to 
explore the rich IR physics of the dual gauge theory. In this paper we focus on the static quark- 
antiquark potential, the screening of chromoelectric charges induced by the dynamical flavors, the 
flux tube breaking and the mass spectrum of the first mesonic excitations. Moreover, we discuss 
the occurrence of quantum phase transitions in the connected part of the static quark-antiquark 
potential. Depending on the ratio of certain parameters, like the flavor mass, with respect to some 
critical values, we find a discontinuous (first order) or smooth transition from a Coulomb- like to a 
linear phase. We evaluate the related critical exponents finding that they take classical mean-field 
values and argue that this is a universal feature of analogous first order transitions occurring in the 
static potential for planar gauge theories having a dual supergravity description. 
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1 Introduction 



The quarks in QCD can be distinguished between hght {u, d, s) and heavy (c, b, t), depending 
on their mass being smaller or (much) larger than the scale A/^j dynamically generated via 
dimensional transmutation. The main vacuum polarization effects related to the ffavors are 
due to the three light quarks, while the others, which can be mainly considered as "probes" 
of the theory, can be neglected in the path integral with a good approximation. 

Nowadays, the only systematic, first-principles-based, non-perturbative approach to low 
energy QCD is provided by the study of extrapolations to the continuum of numerical simu- 
lations of the theory on a Euclidean space-time lattice of finite volume [1] . On the lattice it 
is practically extremely hard to account for the vacuum polarization effects due to the light 
quarks and so many results are obtained using an approximation where they are treated 
as probes. This is the so-called "quenched" approximation. There are certainly indications 
about the possibility of overcoming this limit in the future and some partial results towards 
"unquenching" the lattice are known (see for example [2]). 

In the meantime we have to take in mind that most of what we know from the lattice 
is about quenched QCD. This theory is strictly confining, like pure Yang-Mills, while the 
same is not true for real QCD: due to the presence of the dynamical (light) quarks, an 
external quark-antiquark pair QQ will not experience an indefinitely linear potential at large 
distances. Instead, since a qq pair of dynamical flavors can be popped out of the vacuum, 
the initial QQ state will decay into a pair of heavy-light mesons Qq + qQ for distances larger 
than a certain "screening length" (at which the flux tube becomes sufficiently energetic for 
the decay to happen). 

Remarkably, the string/gauge theory correspondence, which aims to be a complementary 
approach towards explaining the non-perturbative dynamics of QCD, offers a set of simple 
tools to analyze particular unquenched gauge theories in certain regimes. This is certainly 
not an ideal setting at the moment, since the present computational methods allow us to 
give a string dual description to theories which are at most (supersymmetric) extensions 
of planar {N^ — > oo) QCD, where the phenomenologically interesting sector is coupled to 
spurious matter or has some higher dimensional UV completion. Nevertheless, the string 
approach can give valuable insights into many features of strongly coupled gauge theories, 
in regimes inaccessible by other methods. 

In this paper we present an exact (and mostly analytic) fully backreacted D3-D7 super- 
gravity solution dual to the 4d Klebanov-Strassler (KS) cascading gauge theory [3], coupled 
to a large number of dynamical fundamental matter fields which are massive and non-chiral 
(relevant related studies are in Refs. [l]-[23]). The unfiavored theory is known to be confining 
and to share many notable properties with 4d A/" = 1 SYM. The flavored model flows in the 
IR to a SQCD-like theory. The dual supergravity solution allows one to consider, without 
limitations, either light (m < A//?) or heavy (m > A//?) dynamical flavors. Remarkably, 
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the solution is always regular at the origin of the radial coordinates and as such it is a very 
promising tool to explore the IR dynamics of the dual gauge theory. 

With this aim, we begin to explore some relevant physical observables, concentrating 
in particular on their dependence on the flavor parameters. We first focus on the static 
quark- ant iquark potential and study its behavior as a function of the number Nf and the 
masses of the flavors, keeping fixed different relevant field theory scales. We show how 
the solution accounts for the string breaking and screening effects due to the dynamical 
flavors. In particular it enables us to provide a first qualitative study of the behavior of 
the screening lengths as functions of the flavor parameters. We also start up an analysis 
of the mesonic spectrum, by considering fluctuations of the worldvolume gauge field of a 
K-symmetric D7-brane probe corresponding to massless flavors. We show how it is possible 
to study quantitatively how the spectrum varies w.r.t. to the number of sea flavors and their 
masses. Previous studies along these lines in models with flavor backreaction can be found 
in [3 El mm [221 [21]. 

Finally, we analyze the interesting phenomenon of the occurrence of first order quantum 
phase transitions in the "connected" parl0 of the static quark-antiquark potential. These 
transitions are discontinuous changes in the slope of the potential and actually occur in 
several models with and without flavors when at least two separate physical scales are present 
[25| [261 [271 [281 (HI (HI [211 [21] . For the case at hand the potential passes discontinuously from 
a Coulomb-like to a linear behavior. The transition disappears, i.e. the potential becomes 
smooth, for flavor masses above a certain critical value, or by varying some other mass 
parameters in the theory. In the present setup it is possible to evaluate the related critical 
exponents and to show that they take the classical mean-field values. We argue, by means of 
catastrophe theory, that this is a universal feature of every first order transition discovered, 
in different (also unflavored) models, in the static potential holographically evaluated by 
means of string theory. 

1.1 Techniques and structure of the paper 

The addition of fundamental matter to the KS model is realized, following the general sug- 
gestion of [29], by means of space-time filling D7-branes which are holomorphically embedded 
(so as to preserve J\f = 1 supersymmetry) and wrapped on non-compact submanifolds of the 
transverse space [30|. Just as in lattice gauge theory, the task of adding flavors in the stringy 
setup becomes computationally simpler if the vacuum polarization effects due to the funda- 
mental matter fields are neglected. The quenched approximation is realized by neglecting 
the backreaction of the flavor D-branes on the background and treating them as external 

^In the dual picture this refers to the macroscopic open string describing the metastable QQ state when 
the mixing with the final heavy-light mesons Qq + qQ, i.e. the "disconnected" string configurations, is 
artificially turned off. 
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probes. 

In order to go beyond the quenched approximation we use a simple technique which 
was introduced in [6l |7j: we homogeneously smear the flavor branes along their transverse 
directions. This operation is sensible only if the number Nf of such branes is very large, as 
happens in the Veneziano regime (where Nf,Nc — > oo with Nf/Nc fixed). In the smeared 
setups the flavor symmetry group is generically broken to a product of abelian ones, but this 
limitation does not spoil many features of physical interest in the related models. 

When the flavors are massless the related smeared branes reach the origin of the transverse 
space. At this point the flavor symmetry is generically enhanced since the flavor branes 
overlap. In the known cases this is accompanied by the presence of a (good) singularity 
in the dual string solutions. This kind of singularity can be avoided in the massive case, 
where the smeared branes generically extend up to a certain finite distance from the origin 
along the radial direction and there is no special point where the flavor symmetry is fully 
enhanced. For this reason it is extremely interesting to focus on smeared-flavor-brane setups 
where the dual gauge theories are coupled to massive dynamical flavors. 

Using the smearing technique, the addition of massless dynamical flavors to the KS gauge 
theory, giving a solution with a (good) singularity in the IR, was considered in [H]. An 
approximate solution corresponding to the inclusion of massive dynamical flavors to the 
Chamseddine-Volkov-Maldacena-Nunez confining theory [31], was given in [17] and was reg- 
ular in the IR by construction. In order to build up an exact solution accounting for massive 
flavors in the KS case, we have to evaluate the density distribution of the smeared flavor 
D7-branes. This depends on a function of the radial variable, A^/(r), which accounts for the 
effective number of flavor degrees of freedom at a given energy scale. Remarkably, despite the 
complicated setting, we are able to find the explicit expression of A^j(r). The precise knowl- 
edge of this function enables us to write up the fully backreacted D3-D7 solution. Previous 
calculations of the analogous function and the derivation of the corresponding backgrounds 
in the singular conifold case were performed in [THJ [22] . 

This paper is organized as follows. In section [2] we derive the function iVj-(r), the main 
formulas being (12.141) . (12.151) . (12.161) . In section [3] we calculate the full supergravity solution 
dual to the KS model with dynamical massive flavors (the solution can be found in section 
13. 3p . In section m we study the static quark- ant iquark potential and the screening lengths. In 
section [5] we analyze the first order quantum phase transitions in the static quark-antiquark 
potential. In section [6] we focus on some mesonic mass spectra. We end up in section [7| with 
some concluding remarks and a sketch of possible future research lines. The paper includes 
various appendices where many details of the calculations and validity checks are provided. 
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2 Massive non-chiral flavors and smeared D7-branes 
on the deformed conifold 



The deformed conifold is a regular six dimensional non compact manifold defined by the 
equation zi — Z3 z^ = in C^. When the complex deformation parameter e is turned off, 
it reduces to the singular conifold, which is invariant under complex rescaling of the Zi, has 
SU{2) X SU{2) X U{1) isometry and S"^ x topology. The deformation parameter breaks 
the scale invariance, produces a blown- up at the apex of the conifold and breaks the U{1) 
isometry to Z2. 

The low energy dynamics of regular and M fractional D3-branes on the deformed 
conifold is described by a cascading A/" = 1 4d gauge theory with gauge group SU{N + 
M) X SU{N) and bifundamental matter fields A, B transforming as SU{2) x SU{2) doublets 
and interacting with a quartic superpotential Wkw = AiBkAjBi. The KS solution [3] 
is relevant for the N = nM case, where n is an integer. The related theory develops a 
Seiberg duality cascade which stops after n — 1 steps when the gauge group is reduced to 
SU (2M) X SU (M) . The regular KS solution precisely accounts for the physics of an A B- 
symmetric point in the baryonic branch of the latter theory, which exhibits confinement and 
U{1)r Zf2Ar Z2 breaking due to the formation of a gluino condensate (AA) ~ A^^. The 
complex parameter e is the geometric counterpart of this condensate. 

Let us consider the addition of fundamental degrees of freedom to the theory. This can be 
realized by means of suitably chosen D7-branes. A relevant example is given by D7-branes 
wrapping the holomorphic 4-cycle defined by an equation of the form 

zi-Z2 = 2il. (2.1) 

It was shown in [30j that this embedding is K-symmetric and hence preserves the four super- 
charges of the deformed conifold theory. 

A D7-brane wrapping the 4-cycle defined above is conjectured to add a massless (if fi = 0) 
or massive (anti) fundamental flavor to a node of the KS model. The resulting gauge theory 
is "non-chiral" because the flavor mass terms do not break the classical flavor symmetry 
of the massless theory. The related perturbative superpotential is, just as in the singular 
conifold case [32] , 

W = Wkw + hi qM^Bi - As^s)?! + q2{BiAi - 52^2)^2 + ki {q^qif + m (g^g^) , (2.2) 

where we have considered a Z2-invariant setup with two stacks of D7-branes adding the 
same number of fundamental degrees of freedom (with same masses) to both nodes. Thus 
the complex mass parameter m in ly is mapped to the geometrical parameter fi. 

In the singular conifold case the massive embedding zi — Z2 = 2fi explicitly breaks the 
scale invariance and the U{1) isometry of the background geometry. This is related to the 
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explicit breaking of conformal invariance and U{l)ii symmetry due to the mass terms in the 
dual gauge theory. The embedding equation also breaks part of the non abelian symmetry 
group of the conifold to a diagonal SU (2) subgroup. 



2.1 The D7-brane profile 

The metric of the deformed conifold is usually written as 



where 



and 



Kir) 



(sinh(2r) - 2r)i/3 



9 = 



9 = 



9 = 



9 = 



2V3sinh(r) 

— sin 9i d(fii — cos ip sin 02 d(fi2 + sin ip d92 
dOi — sin if) sin 6*2 d(p2 — cos d92 

— sin 9i d(pi + cos ip sin 02 dip2 — sin ip d02 

V2 

dOi + sill ( ' sill 0-2 (1^2 + ' d92 



(2.3) 
(2.4) 



g^ — dip + cos 01 d(pi + cos 02 d(p2 ■ 



(2.5) 



The range of the angles is ip E [0, in), ipi G [0, 2n), 0i G [0, tt], while r G [0, oo). For r — >• oo 
the metric asymptotes the singular conifold one. In terms of these coordinates the non-chiral 
embedding zi — Z2 — 2fi can be written as 



where 



©1 sinh - — i ©2 cosh - — — , 
2 2 e 



^ . Oi . 02 <Pl + <P2-'ip Ol O2 ^l + ^2+1p 

©1 = sm — sm — cos cos — cos — cos 

2 2 2 2 2 2 

. Ol . 02 . ^l + ^2-1p , Ol 02 . (Pl+(p2 + i^ 

©2 = sm — sm — sm h cos — cos — sm , 

2 2 2 2 2 2 ' 



(2.6) 



(2.7) 



and |©i| < 1. From these equations it follows that the profile of the D7-branes has a non 
trivial radial dependence: the branes extend all along the radial direction up to a minimum 
distance Tmin which depends on the relative phase of the p, and e parameters. 
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If fije is purely imaginary (resp. real), the embedding equations imply 9i = (resp. 
62 = 0) and Tmin = Ta = 2arc cosh( / 1 6 1 ) (resp. Tmin = n = 2arcsinh(|/i|/|e|)). For generic 
phases < Tmin < U- Notice that in the "completely misaligned" case where fi/e is purely 
imaginary, Tmin = - i.e. the D7-brane reaches the tip of the deformed conifold - if \fi\ < \e\. 

Let us now define Tq = Ta as the absolute minimal value of Tmin- Related to this, one can 
introduce a parameter mg, the absolute minimum flavor "constituent mass" j§ defined as the 
energy of an hypothetical straight string stretched along the radial direction from r = to 
Tq. If < |e|, then rriq = 0. As discussed above, the dimensional parameters fi and e can be 
related to the bare flavor mass m (see eq. ( 12. 2p ) and the fundamental scale of the dual 
gauge theory. Thus the relation \fi\ < \e\ can be interpreted as m < Ajr. Though at the 
level of the "constituent mass" we do not see differences between the \fi\ < \e\ and the /t = 
cases, we will see that the non zero bare mass parameter influences the density distribution 
of the flavor branes also when m < A/r. This will thus mark a difference with the m = 
setup. 



2.2 The density distribution for smeared D7-branes 

Acting with an S'0(4) ~ SU{2) x SU{2) rotation on the embedding equation zi — Z2 = 2/i, 
we obtain the generalized embedding equation 

pzi - pz2 + qz3 + qz4, = 2fi = 2|/t|e*^ , (2.8) 

where p, q span a unit 3-sphere 

p = cos-e^ 2 >, q = sm.-e^ 2 (2.9) 

and X e [0,47r), G [0,27r), 6 G [0,7r], {3 G [0,27r]. 

Let us now consider a maximal symmetric smeared distribution of Nj ^ 1 D7-branes 
generally embedded as above. By "maximal" we mean that the distribution will not only be 
invariant under the SU{2) x SU{2) isometry of the deformed conifold, but also under the 
symmetry f/(l)^, under shifts of the ip angle, which is broken (to a Z2 subgroup) by the 
deformed conifold geometry. We will thus homogeneously distribute the D7-branes along 
X, 0, as well as along the phase of the mass term j3. We will instead take the modulus 
\fl\ (hence the modulus of the flavor mass parameter in the dual field theory) to be fixed. 
Smearing along (3 will cause different D7-branes to reach different minimal distances r^m 
from the origin. The whole distribution will end up at the absolute minimal distance Tq. 

The density distribution Q of the smeared D7-branes is given by 

^ = J ae,4.,x,f}{6{fi)S{f2)dhAdf2)ded(j)dxd(3, (2.10) 



^This is how the rriq parameter is usually called in the literature. In the present context we adopt the 
same name with an abuse of language. 
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where we have introduced the properly normahzed density function ag^^^^^fj = A^j sin 6'/327r'^ 
and /i = 0, /2 = are the two real constraints implied by the complex equation (12.81) (see 
appendix Rl). 

The symmetries strongly constrain the form ofQ. As was shown in [H], the only possibility 
for an exact two-form preserving SU{2) x SU{2) x U{1)^ x Z2 is, in the present setup, 



Ntir) Nfir) 

' {smeideiAdipi + sme2de2Ad(p2) 7^rfrA(dV^ + cos6'iciv5i+cos6'2rf(y92) , (2.11) 



Air 



An 



where the the dot means derivative w.r.t. r. The function A^/(t) counts the effective number 
of dynamical flavors at a given energy scale, holographically related to r. It crucially depends 
on the particular kind of smeared embedding. Referring to appendix |A] for the details of the 
non trivial and very instructive calculation and defining 



X = coshr, 



(2.12) 



we find that in the present setup the function Nf{x) is the solution of the first order equation 



dNf{x) _ Nffx^e[x-2fx^ + 1] 
dx ~ (x2 - 1)2 

where Q[y] is the Heaviside step function and 



/i(x)-— e [2/^2 + 1 -x]/2(x) 



(2.13) 



h 
h 



ATiil + x"- 



xJ{l + 2/i2 - a;)(l - 2^j? + x){x'^ - 1) + 2/i^(l + x^) arctan 



'l + x){l + 2ij? -x) 



\ (x-l)(l-2^2 + 



X) 



In the massless case /i = we have Nf^r) = const = Nf Moreover, for x < 2/i^ — 1 we 
have Nf{x) = const = 0. 

Let us now split (12.131) into two regions. In region I, x > 2/i^ + 1 (i.e. t > Tb), we find the 
following simple solution 



n\'\x)=Nj 



1 - 2/i^ 



X 



x^ 



(2.14) 



where the integration constant is fixed by consistency so that Nf {00) = Nf. Notice that in 

e^^), this function asymptotes to the expression found 



the large r limit (with [epe"^ 
in [22], for the flavored version of the singular conifold Klebanov-Witten (KW) model [33] . 

In region II, 2/i^ — 1 < x < 2/i^ + 1 (i.e. Tq < r < r^), we have a complicated expression 
in terms of Elliptic integrals of the first and third kind (jF[a | b] and n[a; 6 | c] respectively) 



N 



Nf 



X 



x^ - 1 



TT 



(Ai(x,/i^) +A2(x,/i2) -A2 (2/1^ + 1, 



(2.15) 
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A^/=l , ij=2 





Nf=\ , yU = 0.5 
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Figure 1: Plots of A^/(t) for two examples with n > 1 and < 1. 



where 



(x + - 1) 



X 



— 1 



■ arctan 



2n' 



■if,' - 1)^ 



arcsm 



1 



arcsm 



4 - 1 



(l + a;)(l + 2/i2 -x) 
\ (x-l)(x + l-2/i2) 

+ 



(2.16) 



In (12.151) we have fixed the integration constant by imposing continuity at x = 2/^^ + 1, i.e. 
Nf\x = 2/i2 + 1) = Nf^\x = 2/i2 + 1). This condition is satisfied since Ai(2/i2 + 1, jj^) = 0. 

If /i > 1 and so Xmin = — 1 > 1, we also have that Nf{x) is vanishing for x < Xmin- 
Let us check continuity at x = 2/i^ — 1. First, notice that 



TT (2/^2 - 1) 



so that the above mentioned continuity condition amounts to having 



1 - 



TT 



0. 



(2.17) 



(2.18) 



A2{2fl^ - 1, /i') - ^2(2/1^ + 1, /i2) 

We have checked numerically that this condition is indeed satisfied and, more precisely, that 



1 



TT 



-26 [1 -/i] 



(2.19) 



which thus vanishes for yU > 1. 

Relevant plots of A^/(t) can be found in figure [H For /i > 1 the shape of iVj(r) resembles 
that of a "smoothed-out" Heaviside step function A^/(t) ~ NfQ[T — Tg]. In fact, as we have 
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anticipated, the function iVj(r) counts the effective number of flavor degrees of freedom at a 
given energy scale. At energies smaller than the flavor mass, the fundamental fields can be 
integrated out and the theory resembles the unflavored one. At higher energies the masses 
can be neglected and the theory looks like the massless-flavored one (for which iV/(T) = Nf). 
In other contexts where the density distribution of the flavor branes is difficult to evaluate, 
the Heaviside step function can be fruitfully used to construct approximate solutions. This 
is what was done in [T7] for the massive-flavored CVMN solution. 

Let us stress that also when /i < 1 and so m < Aj^, rriq = 0, there is a non trivial density 
distribution of branes. This marks a difference with the m = case studied in [11]. Of 
course approximating A^/(t) with an Heaviside step function NfQ[T — r^] in this massless 
case would simply replace our solution with the massless one for every r. 

3 The backreacted KS solution with massive flavors 

In this section we present the new supergravity solutions accounting for the full backreaction 
of fractional and regular color D3-branes, as well as of smeared flavor D7-branes on the 
deformed conifold. The solutions follow from an action which is the sum of the bulk type 
IIB supergravity and the flavor brane actions. Following a suggestion in [3^ the action 
for the D7-branes is just taken as the sum of the Dirac-Born-Infeld (DBI) and Wess-Zumino 
(WZ) terms. This is actually an approximation which is sensible only if the effective coupling 
QsNf is small. This is the case in localized setups if Nf <^ or in the smeared setups (where 
Nf can be of the same order of N^) due to the effective suppression of the coupling by the 
large transverse volume [Ml [19j . 

3.1 The ansatz 

In order to present the ansatz for the full background, let us first introduce the one-forms ai 
and cuj (i = 1, 2, 3) as follows 

(Ji = dOi, (72 = sinOidipi, as = cos 6*1 (iy^i , 

oji = sin ip sin 02 dip2 + cos ip d92 , LJ2 = — cos ip sin 62 dip2 + sin ip d02 , 

oj^ = dip + cos92d(p2 . (3.1) 

Here the angles are the same as in the deformed conifold. 

The Einstein frame metric ansatz has the same warped form as in the massless case [11] 



ds' 



h-^/^r) dx^ dxf' + h^/^T) dsl, 



dsl 



le^^3M(rf^2 ^ ^ - g{r)M + gl) + e'^^^^\l + g{r)M + gl) , 
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where dxl^ denotes the four-dimensional Minkowski metric and Gi = Giij) {i = 1,2,3), 
g = g{T) and h{T) are five unknown radial functions. Quite nicely, the embedding equa- 
tion expressed in terms the "deformed conifold r variable" looks the same in terms of the 
"backreacted ansatz r variable" . See appendix [B] for details. 

As for the dilaton and the forms we will adopt the same ansatz as in [11], modulo the 
substitution of Nf with the function iV/(r) evaluated in the previous section. In units gs = ^ 
we have 



B2 
H3 
Fi 
F3 



dh \r) Adx^ A - ■ ■ A dx'^ + Hodge dual . 



a 



a 



2 



fg^Ag^ + kg^Ag' 



Nfir) 

4:71 



drAifg'Ag' + kg^Ag^) + -(k ~ f) g' A {g' A g' + g' A g') 



a -{g' A [{F + -^f)9' A g' + [1 ~ F + 
+FdTA{g'Ag'' + g^Ag^)}, 



An 



k)g'Ag' 



+ 



(3.2) 



where M is the fractional D3-brane Page charge and / = /(t), k = k{T), F = F{t) are 
functions of the radial coordinate (and where the dot denotes derivative with respect to r). 

Notice that, consistently, dFi = —Q, where fl is the D7-brane density distribution form 
given in eq. (12. lip . This and the other modified Bianchi identities 



dF. 



HiAFi - VLAB2 



dF^ = HsAF^ - -QAB2AB2, 
follow from the WZ term of the smeared D7-brane action (see appendix ICll. 



(3.3) 



3.2 The BPS equations 



The modified Bianchi identity for F3 in (13. 3p is automatically satisfied by the ansatz, while 
that for F5 reduces to a first order differential equation for the warp factor 



he 



2G1+2G2 



[/-(/- k)F + ^^fk] + No ^ 



(3.4) 



where Nq is an integration constant that we will set to zero as in [3] and [TT]. The previous 
equation is the same as the one obtained in [11] with the substitution Nf — > A^j(r). In 
general, the BPS equations following from the bulk fermionic supersymmetric variations and 
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from the Bianchi identities are of exactly the same form as those in [TT], with the only 
substitution of Nf with the function Nf{T). This is also due to the fact that, despite the 
modified Bianchi identities of the forms Fi in the massive setup differing from those in the 
massless case, the fermionic supersymmetric variations only contain the Fi and not the dFi. 

In this way one arrives at the same algebraic constraint for g as in [TT]: glg"^ — 1 + 
g2(Gi-G2)j _ g j^g solutions g = and g^ = 1 — e^^^^~^'^^ correspond to the singular 
(and resolved) conifold and to the deformed conifold respectively. Here we focus only on this 
latter case since we want to have a regular solution at r ^ 0. 

The BPS equations for the 6d metric functions are 



Gi = 
G2 = 
G3 = 

while for the dilaton we have 



18 
1 

— ( 
18 



32G3-G1- 



,2G3— Gi— G2 



1 



2 



e 

2 



Gi— G2 



--e 
9 



2G3-G1 



-G2 _j_ gG2— Gi 



in 
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(3.6) 



Notice that, just as in the massless case, by defining Ai = G1 — G2 we get the simple equation 
Ai + 2 sinh Ai = 0, from which, up to an integration constant (that we fix to zero as in the 
massless case) it follows that 



gGi G2 _ ^ g ^ = cosh T . 

Taking this result into account, for the flux functions we have 



(3.7) 
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supported by the algebraic constraint 



, r 
tanh — 
2 



2F coth r + 



Nj{t) 
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k tanh — 
2 



f coth — 
2. 



(3.9) 



In order to solve the above set of equations we will have to distinguish between the two 
possible cases: 1) /i > 1, i.e. m > Am and rriq 7^ 0, with the running function A^/(r) being 
equal to zero for r < r^; 2) /i < 1, i.e. m < and niq = with the running function A^/(t) 
being non trivial up to r = 0. 
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3.3 The solution 



Let us start by considering the \fl\ > \e\ (i.e. /i > 1) case. Correspondingly, there is a region 
T G [0, Tq] where the effective D7-brane charge is zero. In that region, requiring regularity, 
the solution is just (a slight generalization of) the unflavored KS one. Since A^/(t) = 0, the 
dilaton does not run (see eq. (13. 6p ) and the flux functions are just the KS ones, modulo an 
overall constant 

^ rh ^ sinh T — T 

e<P = e'^^« = constant , F = , 

2 sinh T 

f = e<^«^^^^!^(coshr-l)^e'^«/;,5, 
/ smh T 

T COth 7~ — 1 

k = e'^^« {coshT + 1) = e'^"'kKs- (3.10) 

2 sinh r 

The metric is a warped product of 4d Minkowski and the deformed conifold with deformation 
parameter e = em l \2.3L since 



e 



2«i=e2GnanhV, e^^^W = ^^ey^'ir(r) , e2^3W = ^-^. (3.11) 



The warp factor is given by 
2l«''M2 



h{T) 



8 



sinh ^ 



(3.12) 



where h{0) = ho is an integration constant. In [3], /iq was fixed by imposing /i(r = oo) = 0. 
Since the solution we are considering is only valid up to r = we cannot fix the integration 
constant in the same way. 

The above solution has to be continuously glued to the one obtained in the region r > 
where the effective D7-brane charge is non zero. The function iVj(r) has a very non trivial 
expression in general, so we will have to perform some numerical integration. For the dilaton, 
for example, from eq. (13.61) . it follows that 



Air 



NfiOdC, (3.13) 



where the Tq is a point where e*^ blows up. A simple analytic expression can be obtained in 
the r > arccosh(2/i^ + 1) region where A^/(t) is given by eq. (I2.14p 

e<^M = i^- ^ , \ , , — — , (r > arccosh(2u2 + l)). (3.14) 

Nf (ro - r) + 2/^2(1/ sinh(ro) - 1/ sinh(r)) ' ^ K JJ K J 



Just as in the massless-fiavored KS [TT] or in the flavored KW cases [9l [HI [22] , where tq was 
related to a Landau pole in the dual gauge theories, our solution cannot be continued up to 
infinity: r < Tq. 
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Requiring continuity at we find 



-<i>iR 



1 



(3,15) 



which exphcitly depends on Nf and /i. 

A remarkable feature of the present setup is that also in the effectively flavored region 
most of the solution can be given in an analytic way. In fact after some algebra we find the 
following results for the metric functions 



where 



-uv 



e 3 



sinh r 
coshr 



/C(r) 



/C(r) 



g2G. ^£ke-| coshr ^(r), e^^^ = ^ 



[sinh2r - 2r + ^^(r)]' 



/C(r)2 ' 
(3.16) 



(3.17) 



23 sinhr 

and d{4:7ie~'^ri) / dr = (sinh2r — 2r)A^j(r). The function rj^r) is thus a constant in the 
unfiavored region. By requiring continuity at Tg we fix this constant to zero so that 



r/(r) 



47r Jra 



;sinh2e-20iV/(Oc?e. 



(3.18) 



and 



^uv — ^iR e 



<t>IR/'i 



/C(i 



K( 



(3.19) 



'qj V'<?7 • 

The metric in the r > Tg region is thus a warped product of Minkowski 4d and a slight 
deformation (driven by i]{t)) of the deformed conifold metric 



1 I 



efjy e 3 /C(r) 



^ (rfr^ + (9^) + cosh^ (^) ((^7^)^ + (gr) + 



+ sinh^n((,Y + (/)^ 



(3.20) 



In appendix [D] we will verify that this metric reduces to the one found in [TT] when the 
flavors are massless. The warp factor is given by 



h{r) 



AM^a 

8/3 
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di + hi, 
(3.21) 



where the integration constant hi is constrained by the continuity condition at Tg 

(^coth^- l)(sinh2^-20^ 



hi 



8 
^3 

-IR 



ho-e^ 



sinh ^ 



(3.22) 
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Figure 2: Representative plots of the solution for fixed Hq and: yU = 3, Nj = 1 (thick line): 
H = 3, Nf = 1.5 (dashed); /i = 0.5, Nf = 1 (normal). In the plots it is not visible the region 
near the Landau pole, where the dilaton blows up and h ^ —oo. Notice that K. is almost 
not varying with fi and Nf. 



There is something important to notice here: for large values of r tq, h diverges as 
— (tq — r)~^. Thus, /i(ro) = — oo (this happens also in the massless-flavored KS case). Since 
the metric is only well defined for h > 0, we conclude that there exist some maximal value 
of the radial coordinate Tmax < tq where h vanishes and a singularity appears. This behavior 
at Tmax could be connected to the presence of a duality wall [llj. 

For the flux functions things are simpler. By using the constraint equation (13.91) and 
imposing continuity at Tq, we promptly get 



/ 
k 
F 



T coth r — 1 



(coshr 



2 sinh r 

, T coth T — 1 . , , 

e"^ (coshr + 1 

2sinhr ^ ^ 

sinh T — T 



2 sinh r 



(3.23) 



If the D7-branes reach the tip of the cone, i.e. if \fl\ < \e\ (i.e fi < 1), there is no effectively 
unflavored region, and thus the solution has the non trivial A^/(t) dependence in the whole 
r G [0, Tmax < Tq] region. The 6d part of the warped metric is again a generalized (?7(r) 7^ 0) 
deformed conifold with parameter e. In principle in this case we would not need a condition 
analogous to the first equation in (13.191) . In practice, by continuity with the \fi\ > \e\ solutions 
we will rescale e as in (13.191) . 

Sample plots of the relevant functions in the solution are given in figure [21 In order to 
completely verify that these results are correct one has to check that the BPS equations 
solve the equations of motion of gravity, dilaton and forms. To this aim a general result in 
[35] helps: it was shown in that paper that, for general A/" = 1 backreacted solutions with 
metric having a warped (3, 1) x 6 form, the solutions of the first order equations following 
from the bulk fermionic supersymmetry variations and from the (source modified) Bianchi 
identities, are also solutions of the Einstein equations and of the second order equations of 
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motion for the dilaton and the form fields, provided the sources (smeared or locahzed) are 
supersymmetric (i.e. satisfy M = 1 preserving /t-symmetry conditions). Our setup indeed 
satisfies these conditions. In appendix [E] we give the relevant ingredients to perform the 
consistency check explicitly. See also ^14j for relevant related comments. 

3.4 Validity and regularity of the solution 

We can now study the Ricci scalar R and the square of the Ricci tensor RmnR'^^ in the 
string frame, to check the validity and regularity of the solutions we have found above. As 
a first thing, one can check that since e*^ ~ l/^/^ h ~ M'^/Nf (where for h we used the 
rescaling in (13.191) ). we have that R ~ 1/V he^ ~ Nj/M. Thus, as usual in D3-D7 systems, 
the supergravity solution is reliable in the regime 1 -C Nf -C M. The same observation 
applies also to the massless-fiavored KS solution. 

Both curvature invariants diverge at Tmax, though they stay "small" up to values of r very 
close to Tmax (for example up to r = 9.5 if we fix tq = 10)|f| In order to study the behavior 
near r = we need to know the behaviors of the functions h, /C, e"^ in the case /i < 1 (for 
/i > 1 the IR solution is the KS one, so the the solution is certainly regular at r = 0). Clearly 
everything depends on the behavior of A^/(r), which goes as n/r for some constant nj in the 
limit. With this, one can verify that the relevant functions have expansions of the form 

where the explicit expressions for the coefficients (po, 02, ^2, ^755 ^0, ^2 are not relevant here. 
Using these results it is easy to show that the components of the string frame metric have 
expansions in even powers of r around the origin, exactly as in the unfiavored KS case. More 
precisely, our massive-flavored solutions and the KS ones share the same IR topology. It is 
interesting to outline the differences with the singular massless case [H]: there Nf =const, 
g0 ^ q'Po ^ ^ ~ ■?]4 r^, /C ~ /C0 + /C1 r and the string frame metric expansion also contains 
odd powers of r. 

The asymptotics f l3.24p . when inserted in the expressions for R and RmnR^^^ ■, imply that 
the curvature invariants go to a constant at r = and so the solutions are regular in the IR 
for every > 0. 

4 Quark- ant i quark potential and screening lengths 

In the previous section we have constructed a string dual of a flavored version of the 
Klebanov-Strassler theory with a large number Nf of massive dynamical flavors. In this 

■^Note added: there is actually another particular point < Tmax (but very close to Tmax) where the 
curvature invariants are finite but the holographic a-function has a singularity. One can argue the need to 
place a UV cutoff, smaller than Ta, when computing observables in this background [36j . 
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section we begin to study how the dynamical flavors affect the non-perturbative dynamics 
of the gauge theory. We are going to probe the latter with an external quark-antiquark 
Q, Q pair with "constituent mass" Mq. The idea is to study how the static quark-antiquark 
potential, as well as the screening lengths, depends on the sea quark parameters fi, Nj. 

As we have observed in section 12.21 the dynamical flavors have all the same bare mass m 
(related to n) , in modulus, but the corresponding D7-branes have different minimal distances 
Tmin from the origin, since they correspond to different phases of the mass parameter. The 
whole smeared distribution ends up at the absolute minimal distance Tq which we have 
associated to a minimal "constituent mass" parameter rriq. The mass of the probe quarks 
Q,Q is required to be much larger than niq: Mq ^ rriq. The QQ bound state is dual to 
an open string with the extrema lying on a probe D7-brane embedded in such a way that it 
reaches a minimal distance tq Tq, related in the usual way to the "constituent mass" Mq, 
from the bottom of the space0 Due to the existence of a maximal value of r, Tmax < tq in 
the supergravity solution, Mq cannot be taken to be infinite and tq must lie below Tmax- 

The open string attached to the probe D7-brane bends in the bulk and reaches a minimal 
radial position Tup. The Minkowski separation L between the test quarks, as well as the 
potential V{L) (i.e. the total energy renormalized by subtraction of the static quark masses) 
depend on Tup. For an open string embedding given hj t = f,y = a,T = T{y) where 
y G [—L/2, L/2] is one of the spatial Minkowski directions, one finds [37J 



2 

tip 



2 



GP , rQ 



V(Tti„) = / dr - G dr , (4.1) 

27ra'L^, ^P^-Pl^ Jo ^ ^ 

where P, G are expressed in terms of the string frame metric 

P = y/9^=e'^/'h~'/\ G = v/^^=^e^^e*/2^ (4.2) 

and the Hip" subindex means that the quantity is evaluated at r = Tup. It is not difficult 
to check the relation = p^-r^- This for instance means that 4^ and have the 
same sign. In some cases they can change sign simultaneously at some value of Tup so the 
V{L) plot turns around [25j . 

We can now use the background solutions found in the previous sections to study how 
the external quark interaction depends on the dynamical massive flavors. In doing so, we 
vary one of the physical parameters N f, ^{niq), Mq while keeping the others fixed. In the 
numerical analysis below we have set 2^/^M^a'^ej-^''^ = 1 in the expression for h. The masses 



"^More precisely, in order for the corresponding quarks to be non dynamical, we have to take tq much 
larger that the maximal possible value of Tmin in the smeared distribution. As we have discussed in section 
[2] this maximal value is Tb = 2arcsinh(|/i|/|e|). 
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are measured in units of 6~^/^ej|^/(27ra;') and we have put 2'kol = 1. Finally, remember that 
we are working in units Qs = ^ and so what will be denoted by Nf in the following numerical 
studies has to be read as gsNf in standard units. 

As in [T8], we have studied both the exact solutions and those obtained in the Heaviside 
approximation where iV/(r) is replaced by NfQ{T — Tq). Apart from a mismatch in few 
results for the screening lengths, the overall qualitative agreement of the two solutions is 
general. 

4.1 How to compare physical observables 

The solution of section [3] allows to vary the flavor number Nf and mass /x. In general, there 
is no obvious energy scale or coupling which is expected to stay fixed as these parameters 
are varied. Thus, the comparison of various physical observables will crucially depend on 
the choice of the energy scale or coupling which is kept fixed. 

There are actually several possibilities we can consider, since there are many scales in the 
theory at hand. One can decide to fix, for example, the mass of some state in the spectrum 
of the theory. From a lattice perspective a sensible possibility would be fixing the mass of 
the hghtest (0++) glueball@ In our case, this is computationally very hard: the holographic 
calculation of that mass is extremely complicated also in the unfiavored KS background |38j . 
Besides, the lightest glueball excitation in the KS model is actually a massless pseudo-scalar 
[39]. For these reasons, we will not consider this scenario. Another possibility is to keep 
fixed the mass of the lightest meson. We plan to provide in the near future an accurate 
analysis of the mesonic spectrum (which is partly missing also in the unfiavored case), so we 
defer this possibility for a separate publication. We refer to section [6] for some preliminary 
results on the spectrum. Finally, one could also consider the possibility of taking fixed, say, 
the "effective 't Hooft coupling" Me<^ at some scale, for example fixing either its IR value at 
r = or its UV one at the duality wall r = Tmax- 

In this paper we will consider two alternative possibilities, which also have the advantage of 
being computationally easier to realize. They imply different ways of choosing the integration 
constant h^. This constant sets "the scale" of the glueball and KK masses ~ l/m^;„g, 
as can be deduced from the metricl^ Moreover, it enters in the expression for the IR string 
tension T ~ (e'^'^ /hoY^'^ (where e*^^^ explicitly depends on tq, Nj, ji). We will thus consider 
the following possibilities: 

• Possibility 1: constant glueball scale Hq. Keeping constant means fixing the 
glueball and KK scale. The string tension T will change as we change Nf, /z. 

'"'We are grateful to Biagio Lucini and Massimo D'Elia for their coniments on these issues. 
^Clearly, fixing is not enough to fix precisely the numerical value of the mass of the 0++ glueball, 
which depends also on other features of the background. 



18 




Figure 3: The potential at fixed glueball scale Hq = 10. From left to right: Nf = 1, rriq = 1, 
Mq = 40, 50, 60 (dark to pale); Mq = 50, = 1, Nf = 0.6, 1, 1.4 (thick to thin); Mq = 50, 
Nf = 1, nig = 1,5, 10 (continuous, dashed, dotted). 




Figure 4: The potential at fixed IR string tension. From left to right: Nf = 1, niq = 1, 
Mq = 40, 50, 60 (dark to pale); Mq = 50, ruq = 1, Nf = 0.6, 1, 1.4 (thick to thin); Mq = 50, 
Nf = 1, rriq = 1,5, 10 (continuous, dashed, dotted). 

• Possibility 2: constant string tension T. Keeping constant the string tension T 
means performing the rescaling ho h^e'^'^. The glueball scale will thus change as 
we change Nf, fi. 

Let us notice that the string tension is defined from the large distance behavior of the heavy 
(static) quark-antiquark potential in the theory. Since there are dynamical flavors, the heavy 
quark-antiquark bound state is metastable towards decay into a pair of heavy-light mesons. 
Thus one has to keep in mind that for Possibility 2, the string configuration which defines 
T is metastable. 

In both "Possibilities", the radial value tq at which the dilaton diverges (a kind of UV 
Landau pole in the dual gauge theory picture) is kept fixed. Clearly, other choices (other 
"Possibilities") could include the variation of this scale too. We hope to analyze these 
scenarios in the near future. 

In figures [3] and H] we present various plots of the static potential V{L) as the flavor 
parameters are varied using the two prescriptions described above. The potential, which 
has generically a Cornell-like shape, decreases as Mq is increased. The behavior w.r.t. to 
Nf and the mass parameter, instead, depends on the chosen prescription. For Possibility 1 
V{L) decreases with Nf and increases with ruq (with a crossing of potentials for intermediate 
values of L). For Possibility 2 it behaves in the opposite way. 
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Figure 5: The screening length for constant glueball scale ho. In the first plot Mq = 50 and 
the higher (lower) line corresponds to rUg = 1 (/i^ = 0.125). In the second plot Nf = 1 and 
the higher (lower) line corresponds to niq = 1 {fi'^ = 0.125). In the third and fourth plots 
Mq = 50 and Nf = 1. 

4.2 The screening length 

As we have already observed, due to the presence of the dynamical flavors, the QQ state is 
metastable, since a quark-antiquark dynamical pair q, q can be popped out from the vacuum 
causing the decay QQ Qq + qQ. In our setup, the lightest possible heavy-light mesons 
which can arise from the decay of QQ are nearly massless. This is due to the fact that, since 
the dynamical flavor branes are smeared, some of them will intersect the probe one. The 
corresponding heavy-light meson is holographically a string living at the intersection of the 
two branes and its mass is roughly given, at leading order, by Mq/X [JO], where A 1 is 
the bare 't Hooft coupling of the theory. 

The minimal static quark distance at which a pair of these nearly massless mesons can be 
produced is called "screening length" Lg. This is thus defined as 

V{Ls) = -2Mq + 2^ . (4.3) 

Notice that from the relations (14. ip it follows that ^(0) = —2Mq. A non trivial value for 
Lg can thus be obtained only without neglecting the 1/A suppressed mass of the heavy-light 
mesons. We have used a conventional value A = 100 in our numerical analysis, which is 
intended to give a qualitative picture of how the screening length depends on the flavor 
parameterslll The results are in figures [5] (constant glueball scale) and [6] (constant string 
tension) for the two "Possibilities" discussed in section 14. 1[ For consistency with notations 
in previous studies [171 UHl [22], we have used rriq (determined by Tg) as mass parameter in 
the /i > 1 region, where 2/^^ = cosh (r^) + 1. For < 1, on the other hand, rUg is zero, and 
we use directly as flavor mass parameter. 

As anticipated, it is evident from the plots in figure [5] and [6] that the behavior of the physical 
observables in the theory crucially depends on the choice of the fixed scale. The screening 

^We have neglected for simplicity the produced meson interactions as well as the dependence of A on the 
flavor parameters, ft would be interesting to include these contributions in the calculation. Notice that, 
instead, the first a' corrections to the background will induce subleading corrections to the quark-antiquark 
potential. 
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Figure 6: The screening length for constant string tension T. In the first plot Mq = 50 and 
the higher (lower) line corresponds to = 1 (//^ = 0.125). In the second plot Nf = 1 and 
the higher (lower) line corresponds to = 1 (/x^ = 0.125). In the third and fourth plots 
Mq = 50 and A^^ = 1. 

length Lg is in fact a monotonically increasing function of Nf if the glueball scale is kept 
fixed (fig. El first plot), while it is monotonically decreasing if it is the string tension to be 
kept fixed (fig. [6l first plot). For both Possibilities, the screening length Lg is monotonically 
increasing with the static quark mass Mq (figs. [5],[6l second plots). The same is true for the 
behavior with the dynamical flavor mass n in the regime n < 1 (figs. [5], [6], third plots). But 
the behavior for > 1 is again very different: while in the constant glueball scale case Lg 
has a local maximum in the intermediate regime of dynamical flavor masses (fig. [5], fourth 
plot)J^ in the constant string tension case it continues to be monotonically increasing (fig. 
El fourth plot). 

Naively, one would expect Lg to be a decreasing function of Nf (the more flavors there 
are, the more the screening is effective), and an increasing function of /i and rriq (the more 
the flavors are massive, the less the screening is effective). We see that the constant string 
tension "Possibility 2" realizes these expectations. 

Due to the smearing, the decay defining the screening length studied in this section is 
suppressed as l/Nc. In appendix [F] we study a different critical length, the "string breaking 
length", defined in such a way as to reduce this huge suppression to order Nf/N^. The 
results for the string breaking length are qualitatively similar to the ones for the screening 
length, apart from a different dependence on mq. 

5 The quantum phase transitions 

The presence of a quantum phase transition in the heavy quark potential V{L) calculated in 
string theory, appearing as a discontinuity of the first derivative of V{L), is a very common 
phenomenon [25], [261 [211 [281 (HI UHl [211 [21]. It is tempting to conjecture that its occurrence 
is generic in string duals of theories with at least two mass scales, for example confining 

^This feature is not present in the Heaviside approximation, where the behavior is always increasing. 
A possible explanation is that for "sufficiently large" mq the Heaviside is not a very good approximation 
(Nfir) is smoothly increasing). 
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theories with some additional mass scaleo The flavored KS theory is not an exception to 
this behavior. In this section, we discuss the presence of the quantum phase transition 
in the (deformation of the) unflavored case in the massless case of Ref. [llj and in the 
massive case presented in this paper. We also discuss the universality properties of the phase 
transition in all the theories with a supergravity dual, arguing that the critical exponents 
are always classical in the regime of parameter where the gravity description is reliable. 

5.1 Phase transitions in the (flavored) KS theory 

The KS theory is a conflning theory with a dynamical scale Ajr. The integration constant 
ho in the gravity dual was flxed in |3j to some speciflc value hoxs by requiring that h{oo)=0. 
There are indications that a departure of ho from the value hoxs in this type of theories is 
dual to a non-trivial source of a higher dimensional operator [H]. Thus, values of ho different 
from hoKS introduce a second mass scale in the theory. Variations of this scale with respect 
to Aju can cause the heavy quark potential to develop a quantum phase transition. In fact, 
it can be explicitly verifled that this is what happens when ho is tuned above some critical 
value he which is larger than hoxs- Precisely the same phenomenon happens in the CVMN 
solution [31j, as described in [2T]. The features of the potential are always the same in these 
cases and we are going to describe them in the following. This example just indicates that 
also in the flavored case the phase transition is to be expected. 

In the massless-flavored KS case [H], the backreacted D3-D7 solution is singular at r = 0. 
Just as in the massless-flavored CVMN case [7], the singularity causes the static quark- 
antiquark distance L to have a maximum and the V{L) plot to have a turnaround at L = 
Lmax- In flgure [7] we give a sketch of this behavior for a "large" value of the integration 
constant ho- Crucially, for sufficiently small values of ho-, provided the latter is still larger 
than a certain critical value, there is also a second turnaround and a phase transition appears. 

Though the massive case is not continuously connected with the massless one (since our 
gravity solution is always regular), we can flgure out which form we should expect for V{L) 
due to the previous behavior. For /i > 1, the IR solution is just the unflavored one: for this 
reason we see an unbounded linear behavior for the "connected part" of V{L) at large L. 
This behavior can be glued with the massless one allowing the potential to have a double 
turnaround and an asymptotic linear increase. This is actually what it was found in [17] in 
the massive-flavored CVMN theory (the plots are analogous to the ones in flgure [9]). The 
double turn around is such that there is a quantum phase transition between a Coulomb- 

^Confining theories with supergravity duals have actually already at least two distinct mass scales, which 
can be identified with the confining string tension and the glueball scale. Nevertheless, the ratio of these 
two scales cannot be varied at will while remaining in the regime where the gravity description is reliable 
(the string tension must remain much larger than the glueball scale). Due to this reason, for our present 
discussion these two scales cannot be considered independent. 
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Figure 7: The static potential in the massless case for large h^. The IR singularity causes L 
to have a maximal value. At L = L^ax, V{L) turns around. 



Figure 8: Left: the critical value of the flavor mass /ic for the appearance of the phase 
transition, as a function of the number of flavors Nj in the constant IR string tension case. 
Right: in the same case, the V{L) plot for Nf = 1 and increasing values of /i (from left to 
right), showing the disappearance of the phase transition above /be- 
like behavior for the potential at small L and the linear behavior at large L. In |17] this 
quantum phase transition was recognized to be a first order Van der Waals-like transition, 
occurring for flavor mass parameters smaller than a certain critical value. For larger values 
the Coulomb and the linear phase of V{L) are smoothly connected. At the critical point for 
the disappearance of the phase transition, the latter is of second order. 

The solution presented in this paper is dual to a theory with many scales. Besides Am, 
there is a UV Landau pole-like scale related to tq, the scale /i of the dynamical quarks and 
finally the higher dimensional operator scale set by ho- Thus, there is room for a quite rich 
pattern of appearance of phase transitions in the heavy quark potential. In fact, these phase 
transitions appear if the parameters are varied in suitable ways. In some cases, there are 
actually two distinct phase transitions, for example if we keep constant the IR string tension 
and we vary the mass of the dynamical flavors. The first transition happens for mass smaller 
than a critical value fic which is an increasing function of Nf, as can be seen in figure [8] (see 
also figure [9]). The values of /ic can be smaller or larger than /z = 1, so the latter point is 
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Figure 9: The phase transition in the potential V{L) (left) is implied by the presence of a 
local minimum in the Lijtip) (center), or equivalently Vijtip) (right), plot. 

not at all peculiar from this point of view0 The second phase transition happens for mass 
above a second, larger critical value, i.e. when mq approaches Mq. As can be seen from 
the second plot in figure [H which is equivalent to the pressure- volume graph of the Van der 
Waals system, the presence of a local minimum in the Lijtip) plot implies the existence of 
the phase transition. 

It would be very interesting to give an exhaustive description of the patterns of the phase 
transitions in this complicated theory. Our present aim is just to argue that, being the phase 
transition such a generic phenomenon, it is worthwhile to study its properties. In the next 
section we are going to discuss its universality class. 



5.2 Critical exponents, universality classes and cusp catastrophes 

As we have argued above, the presence of the phase transition in our models is signaled by 
a local minimum in the behavior of the quark separation L as a function of the position 
of the tip of the string describing the Wilson loop. Let us focus on the family of functions 
L{Ttip, ii), where we consider for simplicity just the behavior of the flavored KS system as we 
vary the mass /i of the dynamical flavors (equivalently, we could use ho, for example). Near 
the critical point aX fi = fic we have been able to verify numerically, for a selected choice of 
parametera^^l that the function L(rtjp,/i) is well approximated by 

The behavior of the curve near the critical point allows us to extract the related critical 
exponents (for their definition see, for example, [12]) and to determine the universality class 
of our phase transitions. Along the critical line /i = /^c we find that (L/Lc) — 1 ~ [(Ttip/Tc) — 1]*^ 
with the critical exponent S = 3. At L = Lc, fi < fic^e find [{Tup/ Tc) — 1] ~ [1 — (yu/ /ic)]^ with 

^"^Let us note that in the Heaviside approximation one has the very same resuh, even if of course the phase 
transition happens for fic above /i = 1. 

^^We have used gs = I, Nf — 1,to — 10, ho — 10 in the constant IR string tension case. 
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critical exponent (3 = 1/2. The values we have found for S and /3 are thus the classical ones. 
The same results hold for the transition in the static quark-antiquark potential appearing in 
the unflavored KS model when varying Hq. 

Due to the scaling relations, two critical exponents are sufficient to determine completely 
the universality class of the system, once its dimension is known. The scaling relations read 

a + 2/5 + 7 = 2, 7 = P{S-l), 2 - a = ud , 7 = z/(2 - r^) . (5.2) 

The dimension d enters the scaling relations of the critical exponents related to the correlation 
functions, which we cannot calculate; so, the "effective dimension" of our system is undeter- 
mined. On the other hand, the scaling relations let us conclude that, given 6 = 3,f3 = 1/2, 
the other critical exponents are a = and 7 = 1, as in the Landau (mean field) theory, and 
do not depend on the dimension. 

Considering that both in [25] and in [24] the same classical exponents are found for phase 
transitions in different dimensions and driven by different mechanisms, we are led to conclude 
that there must be a universality at work, which would reflect the fact that we are basically 
analyzing a classical object: a macroscopic string. We would like to argue that this is indeed 
the case. 

The universality of the critical exponents is related to the theory of singularity of families 
of functions, or "Catastrophe theory" , see for example [13] and [S] for a related connection 
in the context of charged black holes. 

The basic theorem of catastrophe theory is due to Whitney and was extended by Thom and 
basically states that the singularities of a generic function can be of only a few typesF^ The 
theorem means that a generic function is "equivalent" in the vicinity of the singularity, 
to some "prototype" singularity. The crucial point in this discussion is that the critical 
exponents appear to be invariant among the class of equivalent functions [13] . 

The prototype singularity which interests us is called the "cusp catastrophe" and it is 
deflned by the function 

f{x, a,b) = + ax + b , (5.3) 

which has its critical point at a = 6 = 00. At the "critical isotherm" (a = 0) the condition 
f = x^ + b = determines the value of the critical exponent 6 = 3 from the deflning relation 
b ^ x^ . Analogously, setting 6 = 0, the condition / = + a = determines another critical 

^^Let us stress that our system is "generic" (that is, without special symmetries or properties) by con- 
struction, since we only have it numerically. But in the case analyzed in Ref. |25j the system is described 
analytically and is indeed generic. 

^•^We are actually dealing with families of functions which also depend on some order parameters ai,i ^ 
1, m (for us, it will be to = 2). Two families of functions F, F : M" x M'" M are said to be equivalent 
if there are diffeomorphisms such that F{x,ai) = 

F{y{x,ai),e{ai)) + g{a^). 

-'^''The equation of state of the Van der Waals system is recovered by the map x = V—1, a = T—1, b = P— 1 
in units of the corresponding critical values of the volume V, temperature T and pressure P. 
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exponent (3=1/2 from the defining relation x ~ a^. The corresponding "potential" of the 
Landau theory is just the integral in x of the function / above, and is called "unfolding" 
in the catastrophe theory language. The cusp can be seen by drawing in 3d the 2d surface 
/ = 0, and projecting it on the (a, 6)-plane: in this plane, two curves meeting with a cusp 
at the critical point (a = 0, 6 = 0) are the images of the local maxima and minima of h{x) 
for different a's. For a detailed description of the cusp catastrophe, see for example [43j. 

Remarkably, our curve L^Tup, fi) around the critical point (eq. (15.11) ) is precisely of the 
form f{x, a,b) = provided we map 

x^T^-l ^ b^ — -l. (5.4) 

The curve L{Tup, /i) passes from being monotonic when /i > /i^ to developing a local minimum 
(and a local maximum) when n < fic, where there is a first order phase transition (the tran- 
sition is second order precisely when fi = fic)- For such a family of functions Thom's theorem 
states that the system is equivalent to the cusp catastrophe (15.31) . and this determines its 
universality class. Similar maps hold for analogous phase transitions in the quark-antiquark 
potential in other models with a gravity dual. Thus, the heavy quark potential in the stringy 
regime is equivalent to the cusp catastrophe, every time that there is a phase transition. The 
cusp catastrophe is basically the Van der Waals system, so it implies that the critical expo- 
nents are the classical ones. As said, the crucial point is that the critical exponents appear 
to be invariant among the class of equivalent functions |43]. 

Thus, we are led to conclude that the phase transition in the heavy quark potential in a 
theory with a supergravity dual is in the universality class of the Van der Waals system, at 
least in the regime of parameters where the gravity description is reliable. 

6 A preliminary analysis of the unquenched mesonic 
spectrum 

In this section we consider a probe D7-brane embedded, in the backreacted solution we have 
found, in a similar way as the dynamical flavor brane sources. Fluctuations on the probe 
brane worldvolume are mapped to mesonic modes of spin J = 0, 1. Our aim is to study how 
the spectrum of these fluctuations, hence the mass spectrum of the related mesons, varies 
with the dynamical flavor parameters Nf,fi. 

In general, the task of computing the spectrum, even in the quenched case, is quite difficult. 
This happens because, even if the background metric is simple, the induced metrics on the 
probes can become quite involved and the differential equations for the various modes can be 
non trivially coupled. For simplicity here we will only focus on the simplest fluctuations in 
the probe D7-brane action, i.e. fluctuations of the worldvolume gauge field. The equations 
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of motion to solve will be of the form da{e '^yJ—gF"'^) = 0, where g stands for the induced 
string frame metric. 

Let us start by rewriting our fully backreacted string frame metric, in a form similar to 



e2 



h-^dxl, + B\T){dT^ + (/l3 + h^f) + 



+h-^ A\t) (hl + hl + hl + hl + 



coshr 
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(6.1) 
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h2 
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7 difj „ , 
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Moreover let us rewrite 



hi = (/ii — (^7) cos5 — (/13 sin7 + /i2 COS7) sin 5 
^2 = (Z?-! — (^7) sin5 + (/i3 sin7 + /i2 C0S7) cos5 
^3 = (/is cos 7 — /12 sin 7) + , 



(6.4) 



where we have basically replaced 62, (p2 by 7 G [0, vr], 5 G [0, 47r). 

The reason why we want to make this very complicated change of variables is that the 
non-chiral embedding we are interested in takes a simple form in these coordinates. In 
particular, for the case of a massless embedding zi — Z2 = 0, the D7 profile can be simply 
given by 7 = 0, 5 =const [5U] . 

Following [30j, let us ignore the dependence of the fluctuations on the angles and expand 
the worldvolume gauge field as 



A{x^, t) = e'^"" [a{T)v^dx^ + ai(r)/ii + a^{T)h^)] . 
From the related equations of motion we find (— fc^ = M^) 

dr (A^ tanh rS^a) + M"^ h A^B"^ tanh r a = 

Oi = 



(6.5) 



/coth| \ 
Or — r^Orai + 
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V h 
' A^ tanh r 
B^h 



1 



drdZ + 



-dr{h'^) 
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fi^cothr 
A^h 



+ A^ tanh r 
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Figure 10: Plots of the meson masses with Nf and rriq for fixed glueball (and KK) scale. 
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Figure 11: Plots of the meson masses with Nf and for fixed IR string tension. 

Notice that if we take (j) = 0, these equations reduce to (4.7), (4.9), (4.10) in [30]!^ 

We can now use the standard shooting technique to get the mesonic spectrum. Here we 
only focus on the a(r) mode. We take tq to be fixed and impose regularity at r 0. We 
cannot impose normalizability at r — > cx) since our solutions are defined only up to r^ax 
where h{Tmax) = 0. The most sensible choice seems to treat this point as the position of an 
infinite wall for the fluctuations and so to impose that they vanish at r^ax- These boundary 
conditions render the spectrum discrete. We have examined it for the two Possibilities 
discussed in section HI i.e. fixed glueball mass scale and fixed IR string tension. The results 
are shown in figure [10] and [11] The figures show that, as for the critical lengths, the 
behavior depends on the prescription we choose to fix some physical scale in the theory. For 
Possibility 1 (resp. Possibility 2) we see that the mesonic masses decrease (resp. increase) 
with the number of dynamical flavors and increase (resp. decrease) with their mass. 

The lattice study in Ref. p], performed at fixed bare coupling, shows a decreasing behavior 
with Nf of the difference between the meson masses of two distinct excitations (the n = 3 
and n = 1 modes in our notation). This behavior corresponds to the one of "Possibility 2", 

i^Modulo a typo in equation (4.10) of [30]. 
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i.e. the constant string tension case. Considering also the results for the screening lengths, 
we are led to conclude that keeping fixed the string tension is the prescription giving the 
most natural, possibly useful results. 



7 Summary and discussion 

The study of backreacted supergravity backgrounds offers a big opportunity to explore the 
effects of unquenched flavors in a holographic setup. This analysis could shed light on new 
non-perturbative phenomena that take place in the planar limit of gauge theories. In this 
paper we have developed this program for the case in which the dual gauge theory is the 
conifold theory of Klebanov and Strassler and the flavors are non-chiral and massive, with 
masses either larger or smaller than the dynamical IR scale A^. This is an interesting setup 
to consider, since it amounts to adding dynamical flavors to a conflning gauge theory, like 
we do with pure Yang-Mills to get QCD. 

In our approach, the dynamics of the backreacted gravity plus branes system is governed by 
an action in which the supergravity flelds of the type IIB theory are coupled to the DBI+WZ 
action of the flavor D7-branes. Thus, the D7-branes act as dynamical sources of the different 
supergravity flelds, whose equations of motion and Bianchi identities are modifled by the 
presence of the flavor branes. In the Veneziano limit the number of flavors Nf is large and we 
can homogeneously smear the flavor branes in their transverse space. This has the effect of 
substituting the (5-functions source terms by continuous distributions. Amazingly, the flavor 
distribution function iV/(r), that results from smearing also along the phases of the flavor 
mass terms, can be obtained in analytic form. The smeared D7-branes correspond to flavors 
having all the same mass m in modulus. When m > Ajpt, the distribution extends up to 
a certain flnite value Tq of the radial coordinate. When the mass is smaller than Am the 
distribution extends up to r = 0. It is important to keep in mind that in this case the flavor 
branes do not all extend up to the origin: some flavor branes end up at a flnite distance 
from r = depending on the shift between the phase of the flavor mass term and that of 
the complex deformation parameter e. Thus, differently from the massless case, the origin is 
not a special point where all the flavor branes overlap. 

Knowing the density distribution of the flavor branes one can write an ansatz for the 
supergravity flelds which is a direct generalization of the one for the massless case adopted 
in [11]. The corresponding BPS equations for the different functions of the ansatz, which 
also solve the equations of motion, can be integrated in closed form and the result is a 
generalization of both the unflavored [3j and massless flavored [11] backgrounds. Our new 
solutions have small curvature, and thus are reliable, if 1 <^ Nf -C M. As in the massless 



flavored case of [11] they have an UV singularity at r = Tmax where the warp factor vanishes, 
and a would-be Landau pole at tq > r^ax where the dilaton diverges. However, remarkably, 
the massive flavored solutions are always regular at r = precisely like the unflavored KS 
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solution and contrary to what happens when the unquenched quarks are massless. 

With the flavored solution at our disposal we have started a holographic analysis of the 
non-perturbative dynamics of the dual planar gauge theory. In particular, we have consid- 
ered the properties of bound states formed by an external heavy quark Q and antiquark 
Q. In the presence of dynamical flavors the QQ states are metastable and can decay into 
heavy-light mesons. To characterize these decays we have numerically computed the corre- 
sponding screening lengths and we have studied their behavior as functions of Nf and of the 
mass parameters, keeping constant either the glueball mass scale or the IR string tension. 
We have also presented some preliminary analysis of the mesonic spectrum, obtained by 
analyzing regular fluctuations of a D7-brane probe (corresponding to massless flavors) in the 
backreacted background. 

From the holographic study of the heavy quark potential, i.e. from a classical analysis 
of a static macroscopic open string on the background, we have uncovered the existence 
of quantum phase transitions, which seem to generically occur whenever there are at least 
two distinct physical scales in the dual gauge theory. Our system has several scales and we 
have found numerical evidence that the "connected" part of the static potential undergoes 
first-order quantum phase transitions between a Coulomb-like and a linear behavior. These 
transitions occur until some parameter (like the flavor mass) reaches a critical value, where 
the transitions become of second order. We have evaluated the corresponding critical ex- 
ponents and we have found that, in all the cases, they are given by the classical mean-field 
values. These results, together with those of refs. [25] and [21], led us to conjecture the 
existence of a universality class for the transitions in the static potential of every planar 
gauge theory having a supergravity dual. 

In this work we have just started to examine the predictions of the Klebanov-Strassler 
background with massive unquenched flavor. The fact that we have an IR regular solution 
with new parameters at our disposal (as compared with the unflavored and massless flavored 
solutions) opens new windows for many potentially interesting future studies. 

One problem of obvious interest, for which we have already presented here some prelim- 
inary results, is the analysis of the mesonic spectrum in the backreacted background. We 
would like to determine the dependence of more generic mass levels on the number of un- 
quenched flavors and on their masses, as well as the influence of the different prescriptions 
to determine the parameters of the solution. 

Another interesting direction of future research is the analysis of the entanglement entropy 
for our backreacted solution, following the proposal of ref. jl5] for theories that have a gravity 
dual. In this holographic formulation the quantum entanglement entropy between a region 
A and its complement is obtained as the minimal area of a surface that approaches the 
boundary of A at the boundary of the bulk manifold. This proposal has been applied in 
ref. [IE] to the study of confining (unflavored) backgrounds. The authors of [16] found that 
the entropy of a large class of confining models, including the Klebanov-Strassler solution. 
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displays a quantum phase transition similar to the confinement /deconfinement transition at 
finite temperature. It would be very interesting to analyze how the inclusion of unquenched 
fiavor modifies this results. In particular, by using this entanglement entropy approach we 
should be able to find a phase structure similar to the one uncovered here from the analysis 
of the Wilson loops. 

Finding a finite temperature version of our supergravity solutions would be of great inter- 
est. This would allow us to study the properties of the dual quark-gluon plasma and of the 
corresponding fiavored black holes and to describe how the phase structure of our SQCD-like 
models varies with the fiavor parameters!^ Due to the absence of supersymmetry, the main 
technical problem to tackle in this case is the fact that one has to deal directly with the 
second-order equations of motion of the gravity plus branes system. The stability of a distri- 
bution of smeared fiavor branes in this non-supersymmetric setup is an issue that should be 
analyzed with care. The relevance of this line of research would also be due to the fact that, 
since lattice gauge theory is intrinsically Euclidean, the string/gauge theory correspondence 
is at present the only available tool to explore real-time dynamical properties of strongly 
coupled QCD-like quark-gluon plasmas. 

Finally, it would be important to uncover the details of the field theory which lies in the 
deep IR of the fiavored cascading model, specially in the small (or zero) mass < yU < 1 
cases. Finding the precise maps between field theory vacua and supergravity solutions will 
surely require further investigation. 

We are working on some of these problems and we intend to report on them in a near 
future. 
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A Finding A^/(t) from the smeared embedding 



In this paper we have focused on backreacted D3-D7 solutions on the deformed conifold, 
the manifold defined by the equation ziZ2 — z^z^ = in C4. The D7-brane embeddings 
we have considered are the "non-chiral" ones discussed in [SD], with representative equation 
Zi — Z2 = 2/i. The general family of embeddings along which we have distributed the D7- 
branes, is obtained through an 5*0(4) rotation of the previous equation 



pzi - pz2 + qzs + qzi = 2fi , 



(A.i: 



where p, q span a unit 3-sphere 



6 j(x+jp) 
p = cos -e^ ^2 ' 



. 9 i( x-'p ) 
q = sm -e ^ 2 \ 
^ 2 



(A.2) 



and X e [0,47r), G [0,27r), 6 G [0,7r]. 

To give more explicit expressions in terms of the standard deformed conifold coordinate 
we can use the known relations 
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In the limit r 00, |e| 
singular conifold ones. 

Making use of these formulas we can rewrite (lA.ip as 
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We have defined 

ei = 
6 = 



+ sin 



6i 62 J. ^1 ^2 

sm — cos — cos tz + cos — sm — cos f 4 
2 2 2 2 



cos- 



cos — cos — sm t\ + sm — sm — sm to 
2 2 2 2 



+ 



+ sin 



. Q\ 6*2 . ^ 

sm — cos — sm ^3 



(V^l + V^2 + + ^ + 
- V92 - - ^ + 



6'l . . 

COS — sm — sm 44 



6 

^4 



+ V?2 - + C + ' 
(V9i - V?2 + ^ - ^ + 



(A.5) 



(A.6) 



satisfying tlie relation ^1 + ^3 = ^2 + ^4- Notice also that the 93's only enter through the 

combinations yji + and ipi + Setting 6* = x = = (giving z\ — = 2fi), the A's 
reduce to the 6's in (12.61) . 

By explicit computation, one finds 
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Notice that the ■j/'-dependence has dropped out of this expression. It can be checked that 
Aj + Al< 1. 

From the modulus squared and phase in flA.4l) . one finds the equations 



/i 

/2 



-2 arctan 



sinh^ -Al + cosh^ -Aj - /i' 
2 2 



with 



/i = |/x|e 



2/3 + 47rn = 
= 0, 

i/3 



(A.8) 



(A.9) 



The minimal value of r reached by one of these embeddings is the one for which the inequality 
+ A2 < 1 is saturated. At that point, using the second equation in (lA.Sp . one finds 
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/i^ and A2 



sinh + ji^. Then, from the first expression in (lA.Sp . 
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we find the following relation among the modulus and phase of the mass term and the 
minimal value Tmin 
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A.l Smearing the embedding 



Let us now consider the maximally symmetric smeared distribution of Nf ^ 1 D7-branes 
generally embedded as above. The related density distribution Q, whose general expression 
is given in eq. (I2.10p . can be rewritten, after the trivial integration over /3, as 

N 



^ = ^1 ^^2) sinOdfi A dhdedxd<P . (A.ll) 



Despite this being a quite non trivial integral, we know that the result is severely constrained 
by the symmetries to take the form given in (12. lip . In that expression VLr-^ results to be 
independent on 6'j, ipi and so we can just pick, say, 6'j = = in the corresponding expression 
deduced from (lA.lip . This way we get 

/i = -X - - ^ + 2 arg(6) + 2 arctan [ ~ '''' "| ^ ^ ] -2(3 + Ann , 

— cos + X + (p) 

/2 = 2[coshr-cos(V' + X + 0)]|ercos2^ -4|/i|2, (A.12) 

and so 

{dfi A df2)r,^\^=o = -2[coshr + cos (x + 0)]|e|^cos^ - , (A.13) 

where we have put ip = ^ since Qri> is not expected to depend on ip. Using this expression 
and integrating over 6 we get 

^ Nf /■ 4|/iRcoshr + cos(x + </>)]^,^r , , ,,,, ,2 .,-.,2\ , 

(A.14) 

This has the right limit to the expression found in the singular conifold case [22] when 
r oo,e 0, e'^|ep fixed. Moreover, it is consistent with the embedding reaching the tip 
of the conifold for sufficiently small mass, also in the case where the mass is non zero: the 
Heaviside gives vanishing contribution for small r only if is sufficiently large (compared 
to |e|). In general, the Heaviside fixes the absolute minimal value of r to be given by 
coshTg = — 1 + (2|/ip/|ep) in agreement with the considerations in the previous sections. 
Let us solve the integral above. First of all, let us redefine 

a: = coshr, a = x + 0, l = X~4'i (A. 15) 

such that Xq = 2/i^ — 1. Notice also that | det J\ = 1/2, where J is the Jacobian associated 
to the angular change of variables above. Notice also that requiring 

47r r27v ]^ rmir rnir 

dxd(l)=- / dad'y, (A. 16) 
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implies mn = 16. Hence, let us choose m = n = 4 i.e. [0,47r] as a range for both a and 7. 
Integration over 7 is trivial. From a comparison of (1A.14I) with as expected by symmetry 
arguments (cfr. eq. (12. lip ) we thus arrive at the following integro-differential equation 

dNf(x) Nfii^ (x + cosa) 0, , 

— 1^ = ^^ ^ / Q\x - cos a -2ll^] da. A.17 

dx 2tx Jo (a; - cos a) 2 v/^^^ ^ ^ ^ 

In the massless case fi = pjjj one simply has iV/(r) = const = Nf. Notice that if x > 2^1"^+!, 
the condition imposed by the Heaviside is valid for every a and the calculation gets simpler 
(we can just erase the Heaviside in the integral). The difficult piece of calculation is in the 
region 2/^^ _ 1 < x < 2^^ + 1. After performing the integral in a we get the first order 
equation (12.131) . Solving this equation we get the effective running number of fiavors A^/(t) 
discussed in section 12.21 

The above results have been cross-checked performing numerically the integration in eq. 
(I2.10p . also considering other components of the density distribution form and no contradic- 
tion with the maximally symmetric expression for f2 (in particular the expected ip indepen- 
dence of its components) was found. 



B The holomorphic embeddings on the "backreacted" 
conifold 

All along the paper we have shown that the (Einstein frame) metric ansatz for the back- 
reacted background has a standard warped (3, 1) x 6 form. The 6d transverse space is a 
"fiavor-deformation" , driven by some functions of the radial coordinate, of the deformed 
conifold metric. If in the standard (deformed) conifold case the ^-symmetric embedding 
equations for the D7-branes can be simply written as holomorphic expressions in the Zi, we 
have to ask what happens in the backreacted case. The answer to this question is simple, 
since it can be shown that our backreacted 6d metric is an S'[/(3) structure metric on the 
deformed conifold. Hence the embedding equations for the D7-branes can be written in the 
same way as in the unfiavored case, modulo, eventually, a difference in the radial coordinate. 
To see which coordinate to choose, let us write the deformed conifold as in [48j 

deiW = ~, (B.l) 

where W = Wia^ + W4 and a* are the Pauli matrices. The radial coordinate u of the conifold 
can be defined as 

u^ = i:i{WW+). (B.2) 
It is related to the standard deformed conifold coordinate tdc", by = \e\^ coshroc- The 
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family of Kiihler metrics on the deformed conifold is given in [18] and reads 

du"^ 1 



"■^6 DC 



u 



-A9 



5-\2 



+ ... 



(B.3) 



where we write only the terms which are relevant for the present analysis. Now, our backre- 
acted 6d metric reads 



(B.4) 



and thus requiring it belongs to the family fIB.Sp we get 



9 



'„,2l 



[K"{l-'-^)u^ + K'u 



4 

[K"{1 



|4 



'„.2l 



W + K'u 



dv? 



These conditions can be fulfilled if 

dr^ = ^ 



T = log 



U + \ U 



+ const 



(B.5) 



(B.6) 



which is nothing more than 



u 



^ = leP coshr . 



(B.7) 

if we choose the (irrelevant) "const" above to be: const = — log|ep. From this we see that 
the explicit expression for the D7 embeddings on the backreacted 6d manifold is exactly the 
same as those on the standard deformed conifold: r = tdc- 



C The smeared D7-brane action 

Let us now see which form is taken by the action for the flavor D7-branes in case we homo- 
geneously smear them. The action is taken as the sum of DBI and WZ terms. Let us start 
with the latter, which are simpler. 



C.l The WZ term 



Assuming that the three-form fluxes have only components along the internal directions, the 
WZ term is 



WZ 



Ms 



Cs + B2 A Cq 



A ^2 A (74 



(c.i: 



where Ais is the eight- dimensional worldvolume of the D7-branes, B2 denotes the pullback to 
A^8 of the NSNS two-form B2 and Cs, Ce and C4 denote the puUbacks of the corresponding 
RR potentials of type IIB supergravity. 
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Let us now define the two-form Q as the Poincare dual of A^s; which for any eight-form 



Aa satisfies 



E 



A, 



QAAf^ 



Mio 



In terms of Q, we can rewrite the WZ action of the D7-branes as 

1 



Qsmeared 
'-'WZ 



T, 



D7 



Ml 



Bo AC. 



B2AB2A C4 



(C.2) 



(C.3) 



It is clear from (jC.31) that Vt determines the RR charge of the D7-branes. Notice that the 
wedge product of any form with VL naturally implements its puUback to the worldvolume 
A^8- Moreover, it is clear from flC.3P that VL acts as a magnetic source for the RR field 
strengths Fi, F3 and F5. Indeed, the equations of motion of Cs, Ce and C4 give rise to the 
modifications of the Bianchi identities for Fi, F3 and F5 as given in section [3?T1 



C.2 The smeared /^-symmetry condition 

Let us rewrite the Einstein frame metric ansatz used in the paper, in terms of another radial 
coordinate r, related to r by 3e~'^^dr = dr 



dsl 



dsl 



h{r) ^ dx\ 3 + h{r) 

Let us now introduce the following tangent space basis 
_ ^-1/4 /i = 0, . . . , 3 , e 



h'/'dr 



+ 



-{cus + asfiCA) 



e2 = /iV4eG,^^ 



9(^1) 



g2 ^ /,l/4gG.(^^ ^ 



0G3 



In this basis, the Killing spinors e of the background can be written as 



(C.5) 



(C.6) 



where a is an angle such that tana = —9e^'^ '^^ and is a constant spinor satisfying the 
following projection conditions 



3,1^23,3 {ta )r] = T] 



(C.7) 
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In (IC.7p we have employed a double spinor notation for 77. Notice that these algebraic 
conditions determine four independent spinors, which implies that our configurations are 
1/8 supersymmetric. 

The Kahler form J of the backreacted deformed conifold is the two-form whose components 
are the fermion bilinears J^,y = e{—ia2) T^j^e, where /i, v are indices along the internal 
directions and e is a Killing spinor normalized as ee = h~i. By using (1C.6P and (lC.7p . the 
Kahler form of the transverse 6d manifold can be written as 

/i-l J = e^^+^^ ( (7^ A (7^ - A ^ - ^ rfr A (C.8) 

Making use of the equation ?>{G'i + G'2) = e'^^^^^~'^'^ , which is a consequence of the BPS 
system (IE. 20 . one can check that h~2j is closed 

d{h-^J) = . (C.9) 

The K-symmetry condition for the localized non-chiral embedding is [30] 

JAB2 = 0. (C.IO) 

The smeared version of this condition is 

nAj AB2 = , (c.ii) 

where Q is the D7-brane density distribution form. One can verify that, indeed, this equation 
is satisfied by our ansatz for B2 and the Kahler form. 

C.3 The DBI term 

For zero worldvolume gauge fields the Einstein frame DBI term for the D7-branes can be 
written, in the localized 

Sdbi = -Tdj E / e-^^-det (Gs + e'l B^) , (C.12) 

where the ^'s are coordinates that parameterize the eight-dimensional worldvolume A^g of 
the D7-branes and Gg denotes the puUback to Aig of the ten-dimensional metric G. For a 
background metric with standard warped form, as the one used in this paper, we can rewrite 
the above expression as 

S^j^j = -Td7 T I d^xd^C h-^ Jdet + e't , (C.13) 
where Q = 1, ■ ■ - 4) are the coordinates of the four-cycle wrapped by the branes. 
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The DBI action for the smeared branes takes a remarkably simple form due to the cali- 
bration condition on the wrapped cycles. This condition allows us to rewrite (see eq. (B.IO) 
in [T¥]) 



/det (G4 + e-2 B2) d'^C = \ {J ^ J - e"* B2 A B2) . (C.14) 
Using this result in (]C.13I) . we get 

V/ d'^xh-Ue^ J AJ - B2AB2) . (C.15) 

2 J Ms ^ ^ 



C _ -^^7 
^DBI — 



Following the standard rule, the smeared version of this action is 

gsnreared ^ f d^^^A Vol(Mi,3) A ( c"^ J A J - B2 A B2) . (C.16) 

D Matching with the massless case 



As we have shown in section [2!2| when the flavor mass vanishes, A^/(t) = Nf. In this case 
the integral defining 77 (r) in fl3.18p can be done explicitly and one can verify that K,{t) is 
given by 

/C(r) = .^^^^ , , (D.l) 
43 (ro - r)3 

where A(r) is the function defined in [H], namely 

2(r - ro)(r - sinh2r) + cosh(2r) - 2ttq - 1 
A(r) = ^- — . (D.2) 



sinhr 

Moreover, if we define the new constant e as 



IGtt 

one can verify that the functions Gi, G2 and reduce in this case to: 

or, 1 -.4 sinh^ t . . x 
4 cosh r 

e^*^^ = 7 coshrA(r) , 



[A(r 

which, indeed, are the values found in [TT] . 



4 



39 



E The second order equations and consistency checks 

We verify below that the second order equations of motion for the dilaton, the graviton and 
the various forms are a consequence of the first-order BPS equations. Let us rewrite them 
using the r variable introduced in appendix O For the warp factor the equation is 



4 



f-{f-k)F + ^fk\ + iVo, (E.l) 

where Nq is the integration constant we have put to zero in the paper. For the other metric 
functions we have 

r^' ^ „G3—Gi—G2 '^„G2— Gi— G3 1 ^„Gi— G2— G3 r, 

C.1 -e -e + -e - U , 



lgG3-Gi-G2 I 3^^„_^,_^„ 3 

6 2 



j|lgG3— Gi— G2 _j_ _gG2— Gi— G3 _gGi— G2— G3 



and, for the dilaton 



G' + -e^-^-G^-^^ - 3e^?2-^?i-^?3 + ^ g-^'^^ = , (E.2) 

^e-^-^^ (E.3) 



The BPS equations for the functions k, f and F of the three-forms are 



k' = 3e^~^-^(F + ^f]'-±^ 
\ 47i' J 1-g 



F' = h-f-'^^ik-f) . (E.4) 

It is also interesting to recall that, in the deformed conifold case, the fibering function g is 
related to the metric functions Gi and G2 as 

= 1- e^iG,-G,) 5) 

The equations written above are a consequence of the requirement of supersymmetry. The 
corresponding Killing spinors have been written in appendix O 

E.l The second order equations for dilaton and forms 

The second order equations of motion for the dilaton and the forms in a setup like the one 
considered in the paper are (see also [HI [H]) 

-^aM(G^^^v^9^0) = e'^F^, + l(e^F| - e'^ H!) - ^ 
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d{e^^*Fi) = -e^ifg A*F3 - A B2 A A B2 A Q , 

d{ e't' * F3) = -//g A F5 + i ^2 A ^2 A ^2 A n , 

o 

d{ * H3) = e"^ Fi A *F-i - F5 A F3 + Vol(Mi,3) A ^2 A , (E.6) 

where Vol(Mi 3) is the Minkowski part of the volume element of A^io (for our ansatz 
Vol(Mi 3) = h~^d'^x). In (1E.6P we have not included the F5 equation of motion since it 
coincides with the F5 Bianchi identity which is solved provided the warp factor satisfies the 
first order equation ( IE. II) . 

It is possible to show that the solutions of the first order BPS equations solve the above 
equations of motion. The following hints can be useful in the proof. 

As for the dilaton e.o.m. notice that 



e 



'^F^ - e-^Hl = , (E.7) 



which is actually a consequence of the fact that, as expected from supersymmetry, the 
complex three-form G3 = F(3) + ie~^H is imaginary self-dual in the internal manifold. 
Moreover, it is evident from Sf^f'^'^'^ that the B2 field does not contribute to the equation 
of motion of the dilaton. 

As for the equations for the forms Fi and F^ notice that in our setup 

^2 A ^2 A ^2 A ^2 = , ^2 A ^2 A ^2 = . (E.8) 

One can check that the equation for Fi in ( IE. 61) is satisfied identically, as in the massless 
case. Moreover, but much less trivially, the equation for F3 and are satisfied too due to 
the first order BPS equations. 

E.2 Einstein equations 

The Einstein equations for our system are 

Rmn - ^GmnR = l{dM<pdN<p - ^GMNdp<pd''<p) + ie2*(F£)4^^ - ^GmnF^) + 

I ^ 1^(5) p{5)PQRS ^^0/qp(3) p(3)PQ 1^ pAm 
+ -^^ MPQRS^ N 12 K'^^MPQ^N " l^'^MNi's) + 

e""^ ( 3if|^PQ H^^^ ~ 2 ^'^^^ ) + '^MN , (E.9) 
where Tmn is the DBI contribution to the energy-momentum tensor, namely 

9 12-2 s:Qsmeared 

Tmn - -j^^ . (E.IO) 
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In order to check the fulfillment of (1E.9P it is essential to calculate the different components 
of Tmn- We will compute them by performing explicitly the derivative of the smeared DBI 
action flC.161) with respect to the metric G'^^^ . Let us first consider the case in which M 
and in ( IE. 101) are indices along the Minkowski directions. In this case, the dependence of 
'S'dbT^'^'^ on G^''^" comes from the Vol(Mi^3) volume form. The corresponding derivative is 
straightforward to compute and the result in fiat components with respect to the basis (IC.SP 
is 



1 



7]^^e'^QA\j AJ - e-"^ B2A B2 



(E.ll) 



where the r^'s are the coordinates of the transverse 6d manifold. By using the explicit 
expressions of fl and J (eq. flC.81) ). as well as our ansatz for B2 (eq. (13. 2p ). one can easily 
compute the wedge product of forms appearing on the right-hand side of (IE. lip . One gets 



ill 



-G1—G2 



+ 



Stt 



-G3 



+ 



/ 1,-2 g-^<Jl-^<J2-(-r3 

327r 



kf 



(E.12) 

Let us now obtain the components of Tmn along the internal manifold. Clearly, the only 
dependence of the right-hand side of (1C.16P on the metric of 6d manifold comes from the 
Kahler form J. Then, if /i, u are coordinate indices along the internal manifold and if Ej^ 
are the coefficients of the corresponding inverse vierbein, one haJ^ 



Tab d^V 



nAE^E^-— 



J AJ 



2/iv^ 



^]Ae^ A V [^] A J , (E.15) 



where a, h are fiat indices and e", are frame one-forms along the 6d directions. In a one- 
form basis in which J has the canonical form J = Ae^ + e^ Ae^ + e^ Ae^ it is straightforward 
to verify that the non-diagonal terms A [J] , for a 7^ 6, vanish. Actually, Tab has the 
following diagonal form 



Tab d^rj 



Q A{J A J)e° Sab 



(E.16) 



where {J A J)e<^ is the part of J A J containing the one-form e". It is clear from fIC.Sp that 
one can construct a basis in which J has the canonical form by rescaling appropriately the 



^^Let A^P^ be an arbitrary p-form which, in the basis of the frame one-forms e"^, can be written as 



(E.13) 



Then, we define ie° [A'^^^J as the foUowing (p — l)-form 



(p-1)! 



(E.14) 
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forms g^, - ■ ■ , and dr. After computing the wedge products of the different components of 
J A J with Q, one gets 

nA{J AJ)r = QA {J AJ)g. 



-G1—G2 



IT 



Nf<fr] 



^] A (J A j)^ 



27r 



l,---4) . (E.17) 



It is now straightforward to obtain the remaining components of T in the basis ( IC.Sj) . namely 



T 
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Stt 



fa, 6= 1,2) 



(E.18) 



It is quite non trivial to show, using the previous ingredients, that the solutions of the BPS 
equations also satisfy the Einstein equations of motion. In this verification, the following 
remarkable identity satisfied by the three-forms of our ansatz 



^ ^ MPQ ^ N_ ^ ^ ^MPQ -"jV 



1 

r 



3 o mn , 



M,N = r, 1,2,1,2,3 



(E.19) 



is quite useful. 



F The string breaking length 

As it was observed in [T71 [18], due to the smearing, a generic decay QQ —>■ Qq + qQ is 
suppressed as l/N^ (1/M in our case). To get a decay rate which is not strongly suppressed 
we can consider the possibility of producing a large number of heavy-light mesons up to some 
reference ones. The latter are, arbitrarily, chosen as the mesons whose "dynamical" quarks 
have the same internal charges as the static ones. The string picture of the lowest energy 
configuration for one of these heavy-light mesons is in terms of a straight string stretching 
from the probe Q brane to a dynamical, parallel, q brane. Choosing the dynamical flavor 
brane as the one having the minimal possible distance r^m = Tq from the origin, we get 
that the energy of the string, and so the mass of the corresponding meson, is thus given by 
Mq — nig. The minimal separation at which a pair of such mesons can be produced is called 
"string breaking length" Lsb, and in the present setup0 it is defined as 

r -2mq , /i > 1 , 

ViL,,) = (F.l) 

I , fi<l. 

^^For previous studies of Lgb in a quenched D3-D7 N = 2 setup, see [49] . 
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Figure 12: The string breaking length for constant glueball scale ho- In the first plot Mq = 50 
and the two almost coincident lines correspond to ruq = 1 and /i^ = 0.125. In the second 
plot Nf = 1 and the two almost coincident lines correspond to = 1 and /i^ = 0.125. In 
the third and fourth plot Mq = 50 and Nf = 1. 




Figure 13: The string breaking length for constant string tension T. In the first plot Mq = 50 
and the higher (lower) line corresponds to = 1 (/i^ = 0.125). In the second plot Nf = 1 
and the higher (lower) line corresponds to rriq = 1 (/i^ = 0.125). In the third and fourth 
plots Mq = 50 and Nf = 1. 

The behavior of Lgi, as a function of the flavor parameters is shown in figures [12] (constant 
glueball scale) and [13] (constant string tension), for the two "Possibilities" discussed in section 

The behavior of the string breaking length Lsb is similar to the behavior of the screening 
length Lg, apart from its dependence on rUq, due to the presence of the latter in the defining 
relation (IF. ID . In fact, the string breaking length Lgh is a monotonically increasing function 
of Nf if the glueball scale is kept fixed (fig. [121 first plot), while it is monotonically de- 
creasing if it is the string tension to be kept fixed (fig. [13], first plot). In both cases, Lgb is 
monotonically increasing with the static quark mass Mq (figs. [121 [131 second plots) and with 
the dynamical fiavor mass fi in the regime fi < 1 (figs. [121 [IS third plots). For /i > 1, Lgb is 
monotonically decreasing with /z in the constant glueball scale case (fig. [121 fourth plot). In 
the constant string tension case, instead, Lgb has a very peculiar behavior, displaying both 
a local minimum for small rriq and a local maximum for intermediate values of rriq (fig. [T3l 
fourth plot)0 Let us stress again that, due to the explicit presence of niq in (IF. II) for /i > 1, 
the behavior of Lgb in the latter regime needs not to be continuously connected to that in 
the fi < 1 regime. 

^^In the Heaviside approximation there is just the local maximum, while the local minimum is absent. 
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